Valley polarized electronic beam splitting in graphene 
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We show how the trigonal warping effect in doped graphene can be used to produce fully valley 
polarized currents. We propose a device that acts both as a beam splitter and a collimator of these 
electronic currents. The result is demonstrated trough an optical analogy using two dimensional 
photonic crystals. 
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Since its recent synthesis [1] graphene is proving to be 
a fruitful system in which ideas coming from varied ar- 
eas of physics successfully merge: soft condensed matter 
[2[, pseudorelativistic quantum electrodynamics and, 
more recently, optics |4] and valleytronics The con- 
cept of valleytronics [5] describes the use of the valley 
degeneracy in graphene as a degree of freedom to car- 
rying information in the same way as the spin does in 
spintronic. In the scheme proposed in JB\ the generation 
and detection of the valley polarized current are achieved 
by means of a quantum point contact in a graphene rib- 
bon. The main difficulty in the experimental realization 
of this approach comes from the necessity of tailoring the 
graphene samples into nanoribbons and quantum point 
contacts with zigzag edges needed to support edge states. 
In this letter we propose an alternative mechanism for 
getting a different response from each valley using the 
effect of the trigonal warping [6] that can be tuned by 
a gate voltage in graphene samples of arbitrary shapes. 
The main idea relies on the fact that at energies where the 
trigonal warping is noticeable, the Fermi surfaces around 
each valley are different and anisotropic. This allows us 
to suggest a device based on a p-n junction where the 
response is valley asymmetric. (Figs. 1 and 2). 

In the tight-binding approximation of the low energy 
band structure for graphene, the dispersion relation in- 
cluding both the Dirac and the TW terms, is: 






FIG. 1: (Color online) The two limiting cases for the depo- 
sition of the gate over the graphene sample. The green area 
denotes the region doped with holes at energies where the 
TW becomes relevant, (a) Interface of the gate parallel to 
the zigzag orientation and (b) the armchair interface, (c) and 
(d) represent the corresponding isofrequency lines showing 
the conservation of the parallel wavevector for the zigzag and 
armchair interface, respectively. We obtain two directions of 
refraction for the zigzag orientation parallel to the interface 
depending on whether the chosen valley is K or K' . For the 
armchair interface, a forward direction of refraction appears 
due to the little dispersion produced by the second valley. 
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where vf is the Fermi velocity, a is the lattice constant 
and the short notation /c^ = /c^ + ky is used. The param- 
eter s = ±1 refers to the two valleys. 

The average current moving trough the sample is 
(Jg) = |-|t5pv^^5,where tg represents the transmission 
coefficient and Vg^g stands for the group velocity in each 
valley. When the TW term is considered, Vg^s has the 
form 
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where F = (3A:^-3/c^, —Gkxky). When TW is considered, 
p does not essentially differ from that obtained when 
only the linear term is considered in the whole p — n 
junction 0, 0] having almost the same shape for both 
valleys. The major role in the different behavior of (Jg) 
for each valley is played by the group velocity v^^g, since 
it explicitly depends on s. 

It is well known that optical and electronic systems 
share many wave-like properties. A p — n junction in 
graphene has an photonic analogy with an optical sys- 
tem composed of two effective media with opposite signs 
of their refractive index. In optical material slabs with 
negative refractive index, Snell's law leads to the phe- 
nomenon of light focussing j^, l9|. A similar behavior 
has been proposed in graphene for electrons If the 




FIG. 2: (Color online) (a) Potential profile of a n — p — n~ 
junction and schematic band structure of one valley, (b) The- 
oretical ray tracing of the beams coming from a point source 
in a junction with the interface parallel to the zigzag orien- 
tation. Ef = 0.05eV, Vi = O.lOeV and V2 = -1.5eV. The 
beam splitting is noticeable, (c) Armchair junction. We now 
obtain a collimation perpendicular to the interface. 



Fermi level is chosen such that the conduction band, (for 
which the group velocity Vg points outwards the energy 
surface), in one side of the junction, is connected to the 
valence band in the other side, (where Vg now points 
inwards the energy surface), the system behaves as the 
electrical counterpart of a Veselago Lens [H . 

In Fig{T]we show a field effect n— p+ device in graphene. 
In the area, the applied gate voltage is such that the 
chemical potential is located at energies where the TW 
effect becomes relevant (according to ARPES data 
this appens for energies of 0.6 — 0.7eV). When the TW 
term is considered, the band is no longer isotropic around 
each K point. We have then to care about the relative 
orientation between the lattice and the deposited gate. 
Due to the symmetry of the honeycomb lattice, all the 
possible orientations lie between the two limiting cases 
shown in Figsma),(b). Figdfc) shows the direction of 
Vg^s determined by conservation of the momentum paral- 
lel to the potential barrier interfaces when the boundary 
is parallel to the zigzag orientation. Due to the s depen- 
dence of v^^s the currents associated to each valley travel 
in different directions leading to two different currents 
polarized in the valley index. FiglTJd) represents the par- 
allel momentum conservation when the gate is oriented 
along the armchair direction. The current associated to 
one of the valleys (valley K) is fully collimated because. 
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FIG. 3: (Color online) (a) Transmission probability for the K 
(blue lines) and K'(red lines) as a function of the output an- 
gle Ot for different negative potentials V2. The parameters 
are: barrier width: d = lOOnm, initial energy: Ep = O.OSeF, 
initial potential: Vi = 0.14eV. (b) Output angle Ot of the 
electronic beam as a function of the incident angle 6i for val- 
leys K (upper blue lines) and K' (lower red lines). (See text). 
Barrier width: lOOnm, initial energy Ef = O.OSeV, first po- 
tential step: Vi = O.lOeV. 



by momentum conservation, the accessible area of the 
energy surface is almost plane. The current belonging to 
the other valley is more dispersed because the accessible 
area has a higher curvature. With these characteristics, 
we can design a device to achieve the realization of a 
beam splitter with an n — p — n~ junction (Fig. 2(a)) 
with the interface parallel to the zigzag orientation. In 
the first zone, we introduce an electron source. In this 
region, the Fermi level Ep is at low energy in the con- 
duction band, where the Fermi surface is circular and 
the linear Dirac regime is valid. This area is separated 
from a zone of low chemical potential (p-region) in the 
valence band, by a barrier with a potential Vi . In the 
interface we have both negative refraction and focusing 
[4]. The last step with negative potential V2 is placed 
immediately before the electron focusing in order to get 
a very small illuminated zone. In this way, we have nar- 
rower exit beams avoiding the formation of caustics which 
would appear if the position of the focus were reached be- 
fore the last interface. Finally, the last region is at high 
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FIG. 4: (Color online) (a) Geometry of the trigonal photonic 
crystal (PC) of air holes and optical analogy of a p-n junction 
in graphene. (b) Brillouin zone and isofrequency lines for the 
two PCs for ni (red circles) and for 712 (blue triangles), (c) 
Band diagram for a trigonal photonic crystal of air holes with 
radius r = O.SSapc in a dielectric matrix with refractive index 
m = 3.04 (red bands) and 77-2 = 2.65 (blue bands). The two 
insets show details of the zone around the K point where we 
observe that for the wavelength A = 4.08apc, ni corresponds 
to the conduction band, whereas n2 is located at the TW 
valence band. 



voltages in the conduction band (n~ -region) i.e. it has 
a strong TW. A polarization of valley occurs, this giving 
rise to a beam splitting of the electronic beams belong- 
ing to different valleys (Fig.2(b),(c)). Fig. 3 shows the 
two beams obtained for a potential V2 = —l.beV (TW 
case) 25.4° apart from each other for maximum of trans- 
mission. The beam dispersion is approximately only 3°. 
In the Dirac regime (lower energies) we observe that the 
splitting of the beams is reduced. 

A different situation occurs if the interface is parallel 
to the armchair orientation. In this configuration, (Fig. 
1(b)) the beam of one valley is strongly collimated and 
waveguided (blue color in Fig.2(b)) due to the flat por- 
tion of the isofrecuency line. The other valley slightly 
disperses the beam. Due to the big difference of cur- 
vatures between the Dirac and the TW isoenergies, the 
allowed wavevectors k are situated in a zone of the TW 
isoline with small curvature involving a low dispersion 
(red color in Fig.2(b)). 

Up to day, the most efficient way of identifying 
graphene samples is by optical contrast The maxi- 
mum of this contrast depends both on the thickness of the 
substrate and on the wavelength of the monochromatic 
radiation used to illuminate the sample, usually in the 
visible range. For this range, the maximum contrast cor- 
responds to a thickness of 300nm and the values reached 
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FIG. 5: (Color online) Maps of the electric field modulus of 
a wave originated at a point source located inside a photonic 
crystal with refraction index m and propagating through an- 
other photonic crystal of dielectric matrix 722. (a) Armchair 
orientation, illustrating the beam splitting of the two differ- 
ent valleys K and K'. (b) Zigzag orientation. We can see a 
collimation of the beams since the curvature in the trigonal 
distorted band is smaller than the one of the circular isoline 
in the proximity of the Dirac point. The angular dispersion 
is low. 



for the chemical potential are of the order of O.SeV. How- 
ever, it is theoretically found [l2[ that another maximum 
in the contrast appears at a thickness of 90 — lOOnm. 
With this range, it is possible to get values of the chemi- 
cal potential up to 0.6eV, just in the threshold where the 
TW starts to be appreciable. The effect discussed here 
can then be effectively observed with minor experimen- 
tal changes in current samples, like increasing the gate 
voltage and decreasing the substrate thickness. 

An optical analogy. Control of beam propagation, col- 
limation and focusing in pho tonic crystals (PCs) is al- 



ready well developed p^, [IJ, [151, [161, 117[, thus one can 
assess the validity of the above analysis for graphene by 
putting forward its analogy with this photonic system. 
In fact, an optical analog y o f the graphene system was 
previously described [18|, ll9[. We can realize it with a 
2-D PC with honeycomb, kagome or trigonal structure 
where the gap between the first and second bands is re- 
placed by a Dirac point in the K point of symmetry. As 
an example, we have used a trigonal lattice, with param- 
eter apc, of air cylinders of radius r = 0.33apc inside a 
dielectric matrix with refractive index n (FigHJa)). We 
notice that, in the case of the trigonal lattice, the Bril- 
louin zone (BZ) (FiglJt^b)) coincides with the real lattice, 
while in the graphene honeycomb lattice the BZ is ro- 
tated 30 degrees with respect to the real lattice, which 
implies that the zigzag and armchair situations described 
in graphene are exchanged in the PC. Now we cannot 
change the light frequency throughout the crystal like we 
change the chemical potential in a transistor of field effect 
in graphene, but we can displace the frequency position 
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of the Dirac point in the band diagram by shghtly vary- 
ing the refractive index of the dielectric matrix (Fig|4](c)). 
In this way, a junction of two PCs with the same trigo- 
nal structure but with slightly different refractive index, 
mimics the field effect transistor in graphene. 

Figl5] shows a finite elements numerical simulation 
of the electric field modulus distribution from a point 
source, emitting TE polarized light (i.e. the electric field 
is parallel to the long axis of the air holes) with a wave- 
length A = A.OSapc in a dielectric matrix with refractive 
index ni = 3.04 together with another PC of the same 
structure but a slightly smaller refractive index n2 = 2.65 
for the two limiting cases: armchair ( Fig ISJ a)) and zigzag 
(Figl5jb)). The creation of a beam splitter in the latter 
PC for the armchair situation is illustrated by the fact 
that the two isofrequency curves in the K and K' points 
are triangles pointing in opposite directions from each 
other (right and left). Conversely, if we put the zigzag 
structure parallel to the interface, we obtain a collimation 
of the beams. 

To summarize, we have proposed an efficient way to get 
valley polarized beams using the trigonal warping effect 
in the graphene bands. We have shown a possible design 
of a beam splitter and a collimator device. This forma- 
tion of valley polarized beams, is a key ingredient for the 
use of this new polarization in graphene as an effective in- 
formation carrying degree of freedom. Besides, the trans- 
port properties of graphene can constitute a useful and 
novel device in future nanoelectronics. Finally, we have 
checked the validity of these phenomena, by means of an 
optical analogy with two dimensional photonic crystals 
and proposed a theory of isofrequency lines as a useful 
tool to simulate these kind of processes in graphene sys- 
tems. 
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